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ABSTRACT 

A new turbulence model, based upon dynamic and thermal tur- 
bulent time scale transport equations, is developed and applied to 
homogeneous shear flows with constant velocity and temperature 
gradients. The new model comprises transport equations for k, the 
turbulent kinetic energy; x, the dynamic time scale; k$, the fluctu- 
ating temperature variance; and Xq, the thermal time scale. It offers 
conceptually parallel modelling of the dynamic and thermal turbu- 
lence at the two-equation level, and eliminates the customary pre- 
scription of an empirical turbulent Prandtl number, Pr t , thus 
permitting a more generalized prediction capability for turbulent 
heat transfer in complex flows and geometries. The new model 
also incorporates constitutive relations, based upon invariant the- 
ory, that allow the effects of non-equilibrium to modify the pri- 
mary coefficients for the turbulent shear stress and heat flux. 

Predictions of the new model, along with those from two other 
similar models, are compared with experimental data for decaying 
homogeneous dynamic and thermal turbulence, homogeneous tur- 
bulence with constant temperature gradient, and homogenous tur- 
bulence with constant temperature gradient and constant velocity 
gradient. The new model offers improvement in agreement with 
the data for most cases considered in this work, although it was no 
better than the other models for several cases where all the models 
performed poorly. 
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NOMENCLATURE 

A, B nondimensional parameters 

Cj , c 2 , c 3 , c 4 constants in Z$ equation 

c A , c B constants in c^ equation 

c El , Cg 2 constants in £ equation 

cjl primary coefficient in <v9> equation 

c^ primary coefficient in <uv> equation 

c Tl , c^ constants in x equation 


INTRODUCTION 

Modelling of the temperature field in turbulent flows is rela- 
tively inferior to modelling of the velocity field. The most popular 
methods merely use a constant Prandtl number, Pr t , to infer the 
turbulent thermal diffusivity from a turbulent viscosity computed 
by a two-equation model or an algebraic model. Even in simple 
flows, Pr t is not a universal constant (Reynolds, 1975), and can 
exhibit substantial variation near a wall (Antonia and Kim, 1991; 
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Bagheri, Strataridakis, and White, 1992). Transport models for the 
turbulent heat fluxes have been developed (Shih and Lumley, 
1986; Jones and Musonge, 1988; Lai and So, 1990), but they suf- 
fer from complexities and uncertainties similar to those hindering 
the development and application of transport models for the turbu- 
lent stresses. 

Two-equation models for the turbulent thermal diffusivity have 
been developed (Youssef, Nagano, and Tagawa, 1992; Sommer, 
So, and Zhang, 1993b; Shabbir, 1994) to offer conceptually paral- 
lel modelling of the dynamic and thermal turbulence, and elimina- 
tion of the prescription of Pr t . These approaches permit relatively 
economical, generalized prediction of turbulent temperature fields 
in complex flows and geometries. They employ transport equa- 
tions for k, the turbulent kinetic energy; e, the dissipation rate of k; 
k#, the fluctuating temperature variance; and £$, the dissipation 
rate of kg. In this work, we develop two-equation models that 
employ transport equations for t, the dynamic turbulent time 
scale, and x$, the thermal turbulent time scale, to replace those for 
e and £$. 

The traditional choice of the respective dissipation rates, £ and 
Eg, as the second dependent variables, suffers due to the lack of 
simple, economical boundary conditions. Both dissipation rates 
asymptote to non-zero values at the wall. Exact values can be 
determined from the degenerate form of the k and kg equations at 
the wall, but this involves the computation of the second deriva- 
tives of k and kg, and may cause problems due to excessive stiff- 
ness and additional coupling of the first and second equations in 
each model. 

Our choice of the turbulent time scales as dependent variables 
allows much simpler zero- valued boundary conditions at a solid 
wall. Since x = k / £, and k vanishes at the wall, we can apply the 
same boundary condition for x. Sommer, So, and Zhang (1993a) 
have shown that the assumption of vanishing temperature fluctua- 
tions at the wall is very reasonable, even for specified heat-flux 
boundary conditions. Since x$ = kg / £g, we can apply the same 
assumption for the Tg boundary condition. 

The dynamic k-x model is essentially identical to that developed 
earlier by Speziale, Abid, and Anderson (1992) through direct 
transformation of the k-e equations using the definition x = k / £. 
Our thermal kg-Tg model is developed along similar lines through 
direct transformation of the kg-£g equations using the definition Tg 
= kg / £g. Both models also incorporate constitutive relations, 
based upon invariant theory (Shih and Lumley, 1993), that allow 
the effects of non-equilibrium to modify the primary coefficients 
for the turbulent shear stress and heat flux. 

Although one of the primary benefits of our new model results 
from the simpler boundary conditions to be specified at a wall, it is 
prudent to initially validate any new model by comparison with 
experimental data for simpler homogeneous shear flows. In this 
work, we compare the predictions of our new model with the pre- 
dictions of the models of Youssef, Nagano, and Tagawa (1992), 
designated YNT, and Sommer, So, and Zhang (1993b), designated 
SSZ, with experimental data for decaying homogeneous dynamic 
and thermal turbulence (Warhaft and Lumley, 1978); homoge- 
neous turbulence with constant temperature gradient (Sirivat and 
Warhaft, 1983); and homogeneous turbulence with constant tem- 
perature gradient and constant velocity gradient (Tavoularis and 


Corrsin, 1981). 

In all of these flows, the combination of transverse homogeneity 
of the turbulent correlations and streamwise homogeneity of the 
mean velocity and temperature uncouples the turbulent transport 
equations from the mean transport equations. The homogeneity 
also causes all diffusion terms to vanish, and reduces the turbulent 
transport equations to ordinary differential equations. This greatly 
simplifies the transport equations, and is the primary attraction of 
using homogeneous shear flows for validation of turbulence mod- 
els. 


DEVELOPMENT OF THE DYNAMIC MODEL 

For the simple homogeneous flows considered in this work, the 
standard transport equations for k and e can be written as: 


u l = -< uv >5T e = P ‘- £ 


d£ £ £ 

u d£ = c ^k Pk_CE2 ¥ 


(1) 

(2) 


We note that equation (1) is an exact equation, while equation 
(2) is a generally-accepted modelled equation. 

Using the definition x = k / e, we construct a similar transport 
equation for x (Speziale, Abid, and Anderson, 1992): 


dx x 

u d^ = _c ^k F>k+c « 


(3) 


where c T | = Cgj - 1 = 0.44 and - 1 = 0.8. We note that the 

first term on the right-hand side, representing the transformed pro- 
duction terms, is now negative, given the definition of <uv>, while 
the second term, representing the transformed dissipation terms, is 
now positive. 

Shih, Zhu, and Lumley (1993) derived an algebraic Reynolds 
stress model using invariant theory and realizability constraints. 
For the simple homogeneous flows considered in this work, we 
adopt the linear form of their constitutive equation for the shear 
stress: 


-<uv> = 



(4) 


The primary coefficient, c^, was determined to be a function of 
the nondimensional parameter A: 
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C t* = (1.25+ 1.9 A) 
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The parameter A is proportional to the square root of the ratio of 
production to dissipation, thus allowing to vary from its cus- 


2 


tomary constant value of 0.09, as determined from equilibrium 
and constant-stress assumptions in the log-law region of a bound- 
ary layer. In this region, A = 3.3, leading to = 0.088. 


DEVELOPMENT OF THE THERMAL MODEL 

For the simple homogeneous flows considered in this work, the 
standard transport equations for k# and £$ can be written as: 


*6 d0 

v r* - -< ve >5r £ * ' p « - '* 
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We note that equation (7) is an exact equation, while equation 
(8) is a modelled equation, first suggested by Newman, Launder, 
and Lumley (1981), and later modified by Nagano and Kim 
(1988), to include the second production term. 

Using the definition x e = k$ / E$, we construct a similar transport 
equation for x e : 


dx 0 x e x 0 X 0 

= C xei^ P ke + C te2Y P k + C x63 Y +C tB4 W 

The constants c^i , c^ and c T 04 could be determined 
from c t , c 2 , c 3 , and c 4 in the £9 transport equation. However, since 
the values for these constants in the e e transport equation are not 
as well-tested as those in the standard e transport equation, we 
shall directly determine values for c^, q^, q$3, and c^ by con- 
sidering the behavior of the time scale ratio, R = x e / x, in various 
homogeneous shear flow experiments. 

The new model requires a constitutive equation for the turbulent 
heat flux. The YNT model uses the following relation: 


k 2 

-<v9> = C.-V2R— (12) 

K e dy 

They derived the above form using the nondimensional parame- 
ters A, R, and B, where B is given by: 



We propose an alternative form for B that will be proportional 
to the square root of the ratio of thermal production to thermal dis- 
sipation, similar to A: 


B = 



d0 

dy 



d© 

dy 


(14) 


Using this form, we obtain the following constitutive equation 
for the turbulent heat flux: 


-<V0> = c x kx 6 - (15) 

We note that equation (15) is equivalent to equation (12) if q in 
equation (12) is allowed to incorporate a factor proportional to the 
square root of R (an invariant). 

We show in Figure 1 that the primary coefficient, q, is 
inversely proportional to B 2 for the experimental data of Sirivat 
and Warhaft (1983) for homogeneous turbulence with a constant 
temperature gradient. We also expect that q should have a depen- 
dence on A similar to that proposed for q. Therefore, we propose 
the following form for the primary coefficient to allow it to vary as 
a function of A and B, representing the dynamic and thermal equi- 
librium ratios, respectively: 


k 2 2 d0 

-<ve> = c x7 [2R] 2 - (10) 

while the SSZ model uses: 

k 2 / — d© 

-<v0> = q-j2R— (11) 

K e dy 

In both models, R = x e / x, and q is constant and equal to 0.1. 
Both models are obviously equivalent when R attains its equilib- 
rium value of 0.5, as determined by Beguier, Dekeyser, and Laun- 
der (1978). 

Shih and Lumley (1993) proposed a constitutive equation, 
based upon invariant theory, for the turbulent heat flux. For the 
simple homogeneous flows in this work, its first term reduces to a 
form equivalent to that employed by the SSZ model, although q 
is now allowed to be a function of the invariants: 


From the data of Sirivat and Warhaft (1983), where A = 0, we 
determine c B = 0.77. In order to determine c A , we apply constant- 
stress and constant-flux assumptions, along with dynamic and 
thermal equilibrium, in the velocity and temperature log-law 
region of a boundary layer. In this region, A = 3.3, B = 2.2, and q 
= 0.2, which determines c A = 0.38. 

The constants q^, c^* c te3> ^94 remain to be determined. 
We shall exploit some experimental observations about the behav- 
ior of the time scale ratio, R, using a method similar to that of 
Jones and Musonge (1988). First we consider the decaying homo- 
geneous turbulence experiment of Warhaft and Lumley (1978), 
where R remained constant. We construct a transport equation for 
R from those for x and x$: 
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u— - uT-— e -— — 

dx xdx ^dx 





(17) 


Equating the above to zero determines = 0 and = c t2 = 

0 . 8 . 

Next, we consider the experimental data of Sirivat and Warhaft 
(1983) for homogeneous shear flows with constant temperature 
gradients. Here, the time scale ratio, R, while not constant, 
appeared to approach a constant value downstream. We again con- 
struct a transport equation for R: 


u? = u 

dx 


xdx 


T e dx” 
£ dx 


= c. 


^i^ Pk9+c 
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x03~*2 C t2^2 


(18) 


Equating the above to zero with c^ 3 = c^ determines = 0. 

Finally, we consider the experiment of Tavoularis and Corrsin 
(1981) for a homogeneous shear flow with constant velocity gradi- 
ent and constant temperature gradient Once again, the time scale 
ratio appeared to remain constant through the entire experimental 
domain. We again construct a transport equation for R: 


dR 

= U 
dx 


fl^e x 9 dx 
xdx ^dx 


c x02^ P k + C Te3 


x 2 


+c ^ p *- 


v e 

C x2~ 


x 2 


(19) 


Equating the above to zero with = c x2 determines c te2 = - c xl 
= - 0.44. We can now write equation (9) as: 


dx e T e x e 

U- — — C 1 •- P k. C ^ — — 
dx 1 k k T2 x 


( 20 ) 


We now have a complete thermal turbulence model with only 
two constants which are determined to be equal to their counter- 
parts in the dynamic turbulence model. These constants are related 
to the well-tested standard constants for the e transport equation, 
c el = 1.44, and c^ = 1.8. In contrast, the YNT and SSZ models 
each contain four constants to be determined. The published val- 
ues for all the constants used in the YNT and SSZ models are 
shown in Table 1. 


MODEL VALIDATION 

We consider three classic sets of experimental data (Warhaft 
and Lumley, 1978; Sirivat and Warhaft, 1983; and Tavoularis and 
Corrsin, 1981) to validate our new model for homogeneous shear 
flows. While these are the same data used to evaluate the thermal 
constants in our model, we did not **tune” any of these constants to 


match the data; rather, we used observations about the evolution of 
the variable time scale ratio, as shown by the data, to evaluate 
these constants in terms of the well-established dynamic constants 
used in standard k-e models. 

Various second-moment models (for example: Shih and Lum- 
ley, 1986; Jones and Musonge, 1988; Shikazono and Kasagi, 

1993) have successfully predicted these types of homogeneous 
shear flows. While we cannot hope to compete with a second- 
moment model in predicting these flows, we present our predic- 
tions in an attempt to systematically validate our new model. 

We also compare the predictions of our model against the pre- 
dictions obtained from the YNT and SSZ models, although it 
might be argued that this is an unfair comparison, since models 
developed specifically for wall-bounded shear flows can require 
modification of their constants when used to predict free-shear 
flows or homogeneous shear flows. However, we present these 
comparisons in an attempt to establish a more universally-applica- 
ble model at this level. 

All four transport equations for each model, represented by 
ordinary differentia! equations for homogeneous shear flows, were 
integrated using a standard four-step Runge-Kutta method adapted 
from White (1974). The integration steps were sufficiently small 
to compute results identical to known polynomial and exponential 
solutions of similar equations. Furthermore, computed results for 
the actual model equations using smaller integration steps were 
indistinguishable from those presented herein. 

We first consider the simplest case of decaying homogeneous 
dynamic and thermal turbulence (Warhaft and Lumley, 1978). The 
dynamic turbulence was generated by a uniform biplane grid, 
while the thermal turbulence was generated by a “mandoline” 
array of parallel heated wires.The experiment was conducted at a 
Reynolds number of 10,000 based on a mean velocity of 6.5 m/s 
and a grid spacing, M, of 0.0254 m. Four different initial time 
scale ratios were imposed by varying the downstream location of 
the mandoline, the mandoline wire spacing, and the electrical cur- 
rent applied to the mandoline: Case I, R = 1.04; Case II, R = 0.73; 
Case HI, R = 0.65; and Case IV, R = 0.42. The dynamic turbulence 
field was identical for each of these four cases. All dependent vari- 
ables are normalized by their initial values, while the independent 
variable, x, is normalized by the grid spacing, M. 

The dynamic turbulence results are shown in figure 2. Both k* 
and t* decay at a uniform exponential rate, while x* grows at an 
exponential rate. The predictions of all three models are consistent 
with the experimental data, as expected, since the constant was 
originally evaluated for this type of flow. 

The thermal turbulence results are shown in figure 3, with parts 
(a)-(d) corresponding to the four values for R imposed in the 
experiment. For all four cases, our new model appears to match 
the data more closely than the other models, although they also 
predict an exponential decay for kg* and Eg*, and an exponential 
growth for tg*, as expected. The predictions of the SSZ model 
appear slightly worse than those of the YNT model. 

We next consider homogeneous turbulence with a constant 
mean temperature gradient, but zero mean velocity gradient (Siri- 
vat and Warhaft, 1983; Shih and Lumley, 1986). The dynamic tur- 
bulence was again generated by a uniform biplane grid, while the 
thermal turbulence was generated by either a “mandoline” array of 
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parallel heated wires or a “toaster” array of heated ribbons. The 
mean temperature gradient was maintained by differentially heat- 
ing the individual wires or ribbons. The experiments were con- 
ducted at two Reynolds numbers, based on mean velocity and grid 
spacing; various combinations of heater configurations and 
applied electrical current were then used to generate different val- 
ues of temperature gradient and initial time scale ratio. Table 2 
summarizes the seven test cases used by Shih and Lumley (1986), 
and selected for this work. The dependent variables are normal- 
ized by derived initial scales for each case, while the independent 
variable, x, is again normalized by the grid spacing, M. 

The dynamic turbulence results are shown in figure 4(a) and (b) 
for the high and low Reynolds number cases. The predictions of 
all three models are consistent with the experimental data for both 
cases, and are very similar to those shown in figure 2 , as expected. 

The thermal turbulence results are shown in figures 5 through 
11, corresponding to each of the seven test cases. Parts (a), (b), 
and (c) of each figure show the downstream evolution of k^*, e e *, 
and Xq*; the time scale ratio, R; and the turbulent heat flux, <v9>*, 
respectively. An additional prediction for <v 8 >* from a standard 
k-e model with constant Pr t = 0.9 is shown as a dotted-dashed line 
in part (c) of each figure. 

The predictions of our new model show the best agreement with 
the data for cases 1, 2, and 3, shown in figures 5, 6 , and 7, respec- 
tively. The time scale ratio is adequately predicted as a constant by 
our new model for all cases except case 7, shown in figure 11 , 
which exhibits an initial decay in kg* and R that would be impos- 
sible to model with any gradient-diffusion-based model for a con- 
stant mean temperature gradient. The turbulent heat flux 
predictions exhibit the proper downstream decay in figures 5(c), 
6 (c), and 7(c). The new model predictions for cases 4, 5, and 6 , 
shown in figures 8 , 9 , and 10 , respectively, show poorer agreement 
with the data, especially for the turbulent heat flux in part (c) of 
each figure. 

In all cases, the YNT and SSZ model predictions are in gener- 
ally worse agreement with the experimental data than the predic- 
tions of our new model. The YNT model typically overpredicts 
ke* and <v 8 >*, while the SSZ model typically underpredicts 
them; both models overpredict Tq*. These effects can be traced to 
the nature of the constitutive equations used by each model for 
computing the turbulent heat flux. The YNT model computes 
<v 6 >*, and therefore, production of k^*, as proportional to (2R) 2 , 
while the SSZ model computes it as proportional to (2R)° 5 . As the 
computed value of R departs from its generally-assumed equilib- 
rium value of 0.5, the YNT and SSZ model predictions exhibit 
markedly different behavior. For case 6 , shown in figure 10, R is 
very close to 0.5, and the YNT and SSZ model predictions are 
very similar. For all cases, the prediction of <v9>* using the stan- 
dard k-e model with Pr t = 0.9 is roughly constant, as we would 
expect, and clearly inadequate. 

We finally consider homogeneous turbulence with both mean 
temperature gradient and mean velocity gradient (Tavoularis and 
Corrsin, 1981). The mean velocity gradient (dU/dy = 46.8 s' 1 ) was 
generated by a series of ten parallel channels with individual throt- 
tle screens to vary the mean velocity. The dynamic and thermal 
turbulence was generated by electrically-heated rods placed in 
each channel. The mean temperature gradient (d0/dy = 9.5 K/m) 


was maintained by differential heating of the rods.The dependent 
variables are normalized by derived initial scales, while the inde- 
pendent variable, x, is normalized by the total channel height, h. 

The dynamic turbulence predictions are shown in figure 12. The 
new model predictions for k*, £*, and I*, while not in perfect 
agreement with the data, do show the proper quadratic behavior 
for k* and e*, while x* decays slightly, rather than remaining con- 
stant The new model prediction for <uv>* also exhibits the 
proper behavior, disallowing a slight overprediction. The YNT 
and SSZ model predictions are clearly inferior. 

The thermal turbulence predictions are shown in figure 13. The 
new model predictions for k^*, £ 9 *, and x$* match the experimen- 
tal data very well, as does the prediction for R. The new model 
slightly underpredicts <v 6 >*, although it clearly shows the proper 
trend. The YNT and SSZ model predictions are clearly inferior, 
except for the time scale ratio. The prediction of <v0>* using the 
standard k-£ model with Pr t = 0.9 is also markedly inferior. 

CONCLUSIONS 

A new turbulence model, based upon dynamic and thermal tur- 
bulent time scale transport equations, was developed and applied 
to homogeneous shear flows with constant velocity and tempera- 
ture gradients. The new model comprised transport equations for 
k, the turbulent kinetic energy; x, the dynamic time scale; kg, the 
fluctuating temperature variance; and x^, the thermal time scale. It 
was proposed to offer conceptually parallel modelling of the 
dynamic and thermal turbulence at the two-equation level, and to 
eliminate the customary prescription of an empirical turbulent 
Prandtl number, Pr t , thus permitting a more generalized prediction 
capability for turbulent heat transfer in complex flows and geome- 
tries. 

The dynamic and thermal time scales were chosen as dependent 
variables to eventually exploit their potential for simpler wall 
boundary conditions than those used for the traditional choice of 
the dissipation rates as dependent variables. The new model also 
incorporated constitutive relations, based upon invariant theory, 
that allowed the effects of non-equilibrium in the dynamic and 
thermal turbulence to modify the primary coefficients for the tur- 
bulent stress and heat flux. 

Predictions of the new model, along with those from the YNT 
and SSZ models, were compared with experimental data for 
decaying homogeneous dynamic and thermal turbulence, homoge- 
neous turbulence with constant temperature gradient, and homog- 
enous turbulence with constant temperature gradient and constant 
velocity gradient In general, the new model showed improvement 
in agreement with the data, although in several cases it was no bet- 
ter than the other models. 

The homogeneous shear flows considered in this work repre- 
sented canonical cases that any universal dynamic and thermal tur- 
bulence model should be able to predict Our new model showed 
some ability to offer better predictions than the other models for a 
limited parametric range of constant time scale ratio, constant 
mean temperature gradient, and constant mean velocity gradient 
for these flows. Further work will be required to evaluate the new 
model for the eventual goal of predicting wall-bounded shear 
flows, where the variation of the time scale ratio, dynamic equilib- 


5 


rium, and thermal equilibrium may affect the turbulence differ- 
ently than for these homogeneous shear flows. 
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TABLE 1: SUMMARY OF MODEL CONSTANTS 


Constant 

YNT Model 

SSZ Model 


0.09 

0.096 

c el 

1.45 

1.50 

c e2 

1.90 

1.83 

CX 

0,10 

0.11 

c l 

0.85 

0.90 

c 2 

0.64 

0.72 

c 3 

1.00 

1.10 

c 4 

0.90 

0.80 
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TABLE 2: SUMMARY OF CASES FOR SIRIVAT AND 
WARHAFT EXPERIMENT 


Case 

Re M 

d0/dy 

^■initial 

1 

9545 

10.3 

0.547 

2 

9545 

3.68 

0.610 

3 

9545 

4.48 

0.612 

4 

5150 

1.81 

0.750 

5 

5150 

8.10 

0.772 

6 

5150 

2.24 

0.459 

7 

5150 

1.78 

0.676 


SIRIVAT AND WARHAFT - CASES 1-7 



FIGURE 1: DEPENDENCE OF PRIMARY COEFFICIENT 
FOR TURBULENT HEAT FLUX 



WARHAFT AND LUMLEY - CASE I - R=1 .04 



FIGURE 3(a): THERMAL TURBULENCE PREDICTIONS 
FOR WARHAFT AND LUMLEY CASE I 


WARHAFT AND LUMLEY - CASE II - R=0.73 



FIGURE 3(b): THERMAL TURBULENCE PREDICTIONS 
FOR WARHAFT AND LUMLEY CASE II 



x/M 


FIGURE 2: DYNAMIC TURBULENCE PREDICTIONS 
FOR WARHAFT AND LUMLEY EXPERIMENT 


FIGURE 3(c): THERMAL TURBULENCE PREDICTIONS 
FOR WARHAFT AND LUMLEY CASE III 
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WARHAFT AND LUMLEY - CASE IV - R=0.42 



x/M 


FIGURE 3(d): THERMAL TURBULENCE PREDICTIONS 
FOR WARHAFT AND LUMLEY CASE IV 


SIRIVAT AND WARHAFT - Re M = 9545 
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x/M 


FIGURE 4(a): DYNAMIC TURBULENCE PREDICTIONS 
FOR SIRIVAT AND WARHAFT EXPERIMENT 


SIRIVAT AND WARHAFT - Re M = 5150 



0 60 1 20 180 240 300 

x/M 


FIGURE 4(b): DYNAMIC TURBULENCE PREDICTIONS 
FOR SIRIVAT AND WARHAFT EXPERIMENT 


<QA> 


SIRIVAT AND WARHAFT - CASE 1 
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FIGURE 5(a): THERMAL TURBULENCE PREDICTIONS 
FOR SIRIVAT AND WARHAFT CASE 1 


SIRIVAT AND WARHAFT - CASE 1 
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x/M 


FIGURE 5(b): TIME SCALE RATIO PREDICTIONS 
FOR SIRIVAT AND WARHAFT CASE 1 


SIRIVAT AND WARHAFT - CASE 1 
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x/M 


FIGURE 5(c): TURBULENT HEAT FLUX PREDICTIONS 
FOR SIRIVAT AND WARHAFT CASE 1 








<0A> - 



FIGURE 6(a): THERMAL TURBULENCE PREDICTIONS FIGURE 7(a): THERMAL TURBULENCE PREDICTIONS 

FOR SIRIVAT AND WARHAFT CASE 2 FOR SIRIVAT AND WARHAFT CASE 3 



FIGURE 6(b): TIME SCALE RATIO PREDICTIONS FIGURE 7(b): TIME SCALE RATIO PREDICTIONS 

FOR SIRIVAT AND WARHAFT CASE 2 FOR SIRIVAT AND WARHAFT CASE 3 



FIGURE 6(c): TURBULENT HEAT FLUX PREDICTIONS FIGURE 7(c): TURBULENT HEAT FLUX PREDICTIONS 

FOR SIRIVAT AND WARHAFT CASE 2 FOR SIRIVAT AND WARHAFT CASE 3 
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SIRIVAT AND WARHAFT - CASE 5 


SIRfVAT AND WARHAFT - CASE 4 
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FIGURE 8(a): THERMAL TURBULENCE PREDICTIONS FIGURE 9(a): THERMAL TURBULENCE PREDICTIONS 

FOR SIRIVAT AND WARHAFT CASE 4 FOR SIRIVAT AND WARHAFT CASE 5 


2.0 

SIRIVAT AND WARHAFT - CASE 4 
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SIRIVAT AND WARHAFT - CASE 5 
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FIGURE 8(b): TIME SCALE RATIO PREDICTIONS FIGURE 9(b): TIME SCALE RATIO PREDICTIONS 

FOR SIRIVAT AND WARHAFT CASE 4 FOR SIRIVAT AND WARHAFT CASE 5 


SIRIVAT AND WARHAFT - CASE 4 



SIRIVAT AND WARHAFT - CASE 5 



FIGURE 8(c): TURBULENT HEAT FLUX PREDICTIONS 
FOR SIRIVAT AND WARHAFT CASE 4 
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FIGURE 9(c): TURBULENT HEAT FLUX PREDICTIONS 
FOR SIRIVAT AND WARHAFT CASE 5 









<ve> 


SIRIVAT AND WARHAFT - CASE 6 SIRIVAT AND WARHAFT - CASE 7 



FIGURE 10(a): THERMAL TURBULENCE PREDICTIONS FIGURE 11(a): THERMAL TURBULENCE PREDICTIONS 

FOR SIRIVAT AND WARHAFT CASE 6 FOR SIRIVAT AND WARHAFT CASE 7 



FIGURE 10(b): TIME SCALE RATIO PREDICTIONS FIGURE 11(b): TIME SCALE RATIO PREDICTIONS 

FOR SIRIVAT AND WARHAFT CASE 6 FOR SIRIVAT AND WARHAFT CASE 7 


SIRIVAT AND WARHAFT - CASE 6 SIRIVAT AND WARHAFT - CASE 7 



FIGURE 10(c): TURBULENT HEAT FLUX PREDICTIONS FIGURE 11(c): TURBULENT HEAT FLUX PREDICTIONS 

FOR SIRIVAT AND WARHAFT CASE 6 FOR SIRIVAT AND WARHAFT CASE 7 


11 









<uv> 


TAVOULARIS AND CORRSIN TAVOULARIS AND CORRSIN 



FIGURE 12(a): DYNAMIC TURBULENCE PREDICTIONS FIGURE 13(b): TIME SCALE RATIO PREDICTIONS 

FOR TAVOULARIS AND CORRSIN EXPERIMENT FOR TAVOULARIS AND CORRSIN EXPERIMENT 


TAVOULARIS AND CORRSIN TAVOULARIS AND CORRSIN 



FIGURE 12(b): TURBULENT STRESS PREDICTIONS FIGURE 13(c): TURBULENT HEAT FLUX PREDICTIONS 

FOR TAVOULARIS AND CORRSIN EXPERIMENT FOR TAVOULARIS AND CORRSIN EXPERIMENT 


TAVOULARIS AND CORRSIN 



5 1 0 1 5 20 25 

x/h 

FIGURE 13(a): THERMAL TURBULENCE PREDICTIONS 
FOR TAVOULARIS AND CORRSIN EXPERIMENT 
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